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SUMMARY 

A computer code has been written which utilizes ray-tracing techniques to 
predict the changes In position and geometry of a laser doppler veloclmeter 
probe volume resulting from refraction effects. The code predicts the posi- 
tion change, changes in beam crossing angle, and the amount of uncrossing that 
occur when the beams traverse a region with a changed index of refraction, 
such as a glass window. The code calculates the changes for flat plate, cylin- 
der, general axisymmetric, and general surface windows and is currently opera- 
tional on a VAX 8600 computer system. 


INTRODUCTION 

The laser anemometer has provided the fluid dynamicist with a powerful 
tool for nonintrusively measuring fluid velocities. One of the more common 
types of laser anemometers, the laser fringe anemometer, divides a single 
laser beam into two parallel beams and then focuses them to a point in space 
called a "probe volume" (PV) where the fluid velocity is measured. 

Many applications using this method for measuring fluid velocities require 
the observation of fluids through a window. Whenever light traverses a region 
of a different refractive index, its direction of travel is changed in a man- 
ner described by Snell's law: 


N sin a =» N'sin a 1 


( 1 ) 


where 

N,N' indices of refraction before and after the surface 

a, a' angles between the surface normal and the light beam on either side of 
the surface. 

With the exception of a few very specific window/optics configurations, 
the Implications of this law to the laser probe volume when a window is 
inserted In the laser beam optical path between the focusing lens and the probe 
volume can be summarized as follows: 



(1) The position of the actual probe volume will change. 

(2) The beams will uncross, i.e., no longer lie in the same plane at the 
probe volume location. 

(3) For nonflat plate windows, the crossing angle of the two beams will 
change. 

The first result means that data will not be acquired in the desired loca- 
tion. The second affects the data acquisition rate by distorting the probe 
volume, which only exists where the beams cross. The last impacts the 
accuracy of the measurements taken because it is directly related to the calcu- 
lation of velocity by the following equation: 


V = Xf/(2 sin(e/2)) 


( 2 ) 


where 

X light wavelength 
f fringe crossing frequency 
e beam crossing angle 

In the past, most researchers have assumed that the resulting errors are 
small and can be Ignored. The few attempts to assess the seriousness of win- 
dow refraction effects (refs. 1 to 4) have, unfortunately, been concerned with 
special cases such as flat windows and cylinders with observations in a single 
plane normal to the axis of symmetry. This approach has proved reasonable in 
the past since windows were usually simply shaped and thin. Further, observa- 
tion angles were generally near normal to the window surfaces. 

A more general approach developed at the NASA Lewis Research Center to 
meet the needs resulting from the use of more complex window shapes and pas- 
sage geometries is reported here. 

A ray tracing technique which is not restricted to special cases is used 
to study the changes in probe volume geometry and position due to refraction 
effects caused by both flat and general smooth windows. This technique pre- 
dicts the laser beam paths through any of four different window shapes: 

(1) Flat plate windows (constant thickness) 

(2) Cylindrical windows (constant thickness) 

(3) General axisymmetric windows (constant thickness) 

(4) General smooth windows (variable thickness) 

When describing the window surfaces, the technique calculates the flat 
plate and cylindrical window cases analytical ly. The general axisymmetric 
case is treated analytically in the R-6 plane and uses a cubic spline fit in 



the R-Z plane where the Z axis represents the axis of window symmetry. 

The smooth general window case uses the cubic spline fit in both the X-Y 
and X-Z planes. 

The technique yields the new probe volume position, the new crossing 
angle, and the amount of beam uncrossing (specified as both the number of 
probe volume fringes and the absolute minimum distance between beams). 

Two separate computer programs are used. The first program generates an 
input data set for use by the second program. The second calculates the 
changes in the geometry of the probe volume. The use of a separate program 
allows the creation of permanent input datasets. 

The following parameters may be varied to describe the case to be 
evaluated: 

(1) Indices of refraction on both sides of the window and the window 
Itself 

(2) The window shape 

(3) The assumed position of the probe volume and the actual position of 
the focusing lens relative to the window 

(4) The orientation of the plane which contains the laser beams 

(5) The beam crossing angle 

(6) The laser beam wavelength 

This report contains a description of the problem geometry, the method of 
analysis, complete Information for using these programs, example input and out 
put, and several example cases for general axisymmetric and smooth totally gen 
eral surfaces. 

These FORTRAN 77 programs are currently operational on a DEC VAX8600 com- 
puter. These Interactive programs require a typical core memory of 86 pages 
and have typical CPU run times of less than 1 sec, depending upon window 
geometry. 


GEOMETRY 

A basic understanding of the window and optics geometry and the variables 
used to describe them is necessary to correctly use the technique explained in 
this report. Therefore, we will take a more detailed look at the geometries 
of the surfaces involved. Figure 1 provides Information on the general orien- 
tation of the probe volume and lens for all windows. While most of the infor- 
mation presented in this figure is self explanatory, some points should be 
clarified. 


The X', Y' t V coordinate system shown in figure 1 has its origin at 
the probe volume. The X-Y-Z coordinate system, not shown in figure 1, is 
oriented such that the Z axis is along the axis of symmetry of the window 
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being described, should this window be symmetric. The axes of the X', Y', 

V coordinate system are parallel to their respective X, Y, Z axes. The 
origins of the two coordinate systems are not necessarily coincident. 

The X axis should be oriented roughly so that the optical system bisect- 
ing line parallels the X-Z plane. 

The bisecting line is defined as a line equidistant from both laser beams 
and lying in the plane defined by them. It connects the lens center and the 
original probe volume. 

The original probe volume is the probe volume position that would exist 
if there were no window in the optical path between the focusing lens and the 
probe volume. 

The beam crossing angle, e, is, for all cases, the angle defined as: 


tan e = | A-| x A 2 | /A] • A 2 


(3) 


where 

e beam crossing angle 

A], A 2 incoming laser beam vectors 

The angle 4>, not shown in figure 1, is the angle between the surface nor- 
mal and the bisecting line for the flat plate window. 

For all cases, the beam plane orientation angle, 5, is defined as the 
angle measured counterclockwise in the lens plane, made by the line connecting 
the two beams where they cross the lens plane and a line, Z", where Z" is a 
line lying in a plane formed by the V axis and the beam bisecting line 
(fig. 1). 

Figure 2 shows the geometry for the cylindrical window case and additional 
information for describing the lens and original probe volume position. Notice 
that the axis of symmetry of the window is the Z-axis. 

For all window cases, the lens center and probe volume position are 
described in the R, Z, © coordinate system. Neither the probe volume nor 
the lens center axial position are restricted to a Z equal zero requirement. 

Figure 3 shows the general axisymmetric coordinate system. The window is 
described as a set of R-Z points that specify the window outer surface. The 
window inner surface is specified using a normal thickness. For this window, 
the thickness must be constant. 

Figure 4 indicates the method in which the general smooth window is speci- 
fied. Here both surfaces are specified as an array of X, Y, Z points in 
space. There are no restrictions on window thickness in this case. The num- 
bers in parentheses indicate the array locations of the points. 
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The specific restrictions to the geometries described above are: 

(1) Surfaces must be smooth and smoothly changing. 

(2) For all general axisymmetrlc surfaces, surface slopes In the Z 
direction, DX/DZ, must remain less that 45° and DY/DZ must remain less than 
45° for all general smooth surfaces. 

(3) The Z axis Is defined as the center of rotation for the cylindri- 
cal and general axisymmetrlc surface. 

(4) R1 , the lens radial position, should be greater than Rpv, the probe 
volume radial position, and greater than Ro, the window outer radius and 

Rpv should be less than R1 , the window Inner radius (fig. 2). 


METHOD OF ANALYSIS 

The code uses a simple ray tracing technique for the necessary calcula- 
tions. To accomplish this, the initial ray direction cosines and the points 
In space where the beams cross the lens plane are calculated. 

Then, from the window surface information, the point of Interception of 
each beam with the upper or outer window surface and the surface window normal 
at that point are calculated. Snell's law is applied to each beam as It 
crosses the surface. The same procedure is then applied to find the intercep- 
tion of the beam and the lower or inner window surface. 

With the position and the direction cosines of the two beams at the inside 
of the window, the calculation of minimum crossing distance, crossing angle, 
and the point of nearest approach become relatively straightforward geometry 
problems . 

Of importance is how the beam/window interception point and the surface 
normals at the interception point are calculated. 

The following section discusses the methods used to determine the point 
of beam/window interception and the window outward surface normal at the 
beam-window interception point for each of the four window types this computer 
code solves. 

At the point of intersection of the laser beams and surfaces, the laser 
beams and the surface normals can be described using direction cosines in the 
following way: 


V = (cos a) i + (cos (3) j + (cos y)k (4) 


where 

V incoming laser beam vector 

a, (3, y direction angles between the vector and the X, Y, Z axes 
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for each incoming laser beam and 


N » (cos u)i + (cos 03 + (cos o)k 


(5) 


where 

N outward surface normal 

v,C ,o direction angles between the vector and the X, Y, Z axes 

for each outward surface normal (figs. 5 and 6), where both are defined as unit 
vectors. 

The angle between these two vectors is simply the arc cosine of the dot 
product. 


COS* 1 ^ = V • N a COS a COS u + COS |5 COS C + COS y COS a 


( 6 ) 


Using a simple ray tracing analysis and Snell's law; the direction of the 
ray after crossing the surface can be determined by (ref. 5): 


V' = (n/n ' )V + (b/n 1 )N 


(7) 


where n is the index of refraction on the incoming side of the surface, n 1 
is the index of refraction after crossing the interface, N is defined by 
Eq. 5, and b is: 


b = n 1 cos <t> - n cos <t>' 


( 8 ) 


and <J>' is given by 


cos <t> * = (n'2 _ + n^ cos^ a)/n‘ 


(9) 


Using these equations, a method for describing the window surfaces in 
space, and the window surface normals; the path of each beam can be traced. 

The window surfaces and their surface normals are described differently 
for each of the four types of surfaces. 

The flat window surfaces are described using the equation of a flat 
surface: 


(X - Xl)/cos v = (Y - Yl)/cos C = (Z - Zl)/cos o (10) 
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where 


X, Y,Z surface coordinates 

XI. Y1.Z1 a specified surface point 

v,C,a surface normal direction cosines 

A flat surface can be specified by any three points not In a straight 
line. The cross product of any two unit vectors connecting these points pro- 
vides the surface normal and Its direction cosines. Since the surface Is 
flat, the surface normal has a constant direction. 

The cylindrical and general axlsymmetrlc windows are described using the 
equation: 


r 2 * x 2 + Y 2 (11) 


where 

R constant radius for a cylinder and 

R f(z) for the general axlsymmetric case and f(z) Is specified by a series 

of R-Z points 

X,Y surface coordinates 

The necessary direction cosines for a cylinder are also relatively 
straight forward calculations. The slope of the circle at any point can be 
obtained by taking the derivative of X with respect to Y. This Informa- 
tion, coupled with the fact that the direction cosine with respect to the axis 
of symmetry is zero provides all the necessary information. 

To approximate the unknown R-Z relationship in the general axisymmetric 
case, the program makes use of a cubic spline fit subroutine described In ref- 
erence 6. It solves a tridiagonal matrix given in reference 7. The Inputs 
are a series of R-Z values. It should be noted that the spline subroutine 
requires that the end point condition of the second derivative at either curve 
end point be specified. The Input default gives an end point second deriva- 
tive of one half the second derivative of the adjacent spline point. Care 
must be taken to Insure the Input points accurately describe the surface, that 
the surface is smooth and that the end point second derivatives are not greatly 
In error. Inaccuracies occur In the solution when the surface described has 
slopes approaching 90° and surfaces of this type should be calculated using 
the general surface routine with an appropriate coordinate rotation for the 
geometry. 

Output from this cubic spline fit at the desired Z location includes 
the value of the radius and both the first and second derivatives of radius 
with respect to axial position. To solve for the intersection between the win- 
dow surface and incoming laser beam, the Incoming beam Is first described as a 
function of radius and axial position. An Iterative procedure is then used to 
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find the R-Z interception point of the incoming ray and the surface. The de- 
sired solution accuracy for the Iterative procedure is an input value. 

At the calculated intersection point, the slopes of the general axisymmet- 
ric surface are known: 

dR/dZ and dX/dY 

However, we are interested in the surface normal direction cosines (fig. 5): 
cos v, cos C. and cos a. The last is directly available since: 


tan a = -1 /(dR/dZ) (12) 


Now, recalling that: 


1 = (COS 2 v) + (COS 2 C> + (COS 2 a) (13) 


and that, since dX/dY is known 


tan 5 = dX/dY (14) 


Now, 


cos v = b/s (15) 


and 


s = (a 2 + b 2 + d 2 ) 1/2 (16) 


and realizing that 


a = b tan £ and d = b tan k (17) 


where a, b, c, and d are distances defined in figure 5(b) we can write 


cos v - b/(b 2 + b 2 tan 2 £ + b 2 tan k) 1/2 (18) 


or 
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cos v = 1/(1 + tan 2 £ + tan 2 *) 1 / 2 


(19) 


and 


cos 5 ■ (1 - cos 2 a - cos 2 v) 1 / 2 


( 20 ) 


The general surface solution follows the same logic; the difference being 
that no analytical expressions exist to describe the surface. The surfaces 
are described by an array of X, Y, Z coordinates. These coordinates should 
form a roughly rectangular mesh (fig. 4). The method uses a linear Interpola- 
tion scheme (ref. 6) to determine an approximate ray-surface Intersection 
point. Once a "near solution" point Is found, the method calls the spline fit 
routine In the two orthgonal directions (Y and Z In fig. 6) to calculate a 
more accurate Interception point. The subroutine checks the Interception 
points and the changes In slope using the spline fit In both orthogonal direc- 
tions. Should they be outside the Input tolerances, an error message Is 
printed. However, the calculation continues. 

The surface normal Is calculated In a manner similar to that of the 
general axisymmetrlc case. The differences are that the outputs from the 
cubic spline fits Include: 


dX/dY and dX/dZ 


Thus 


tan £ - 1 /(dX/dY) 


( 21 ) 


and 


tan K - 1 / (dX/dZ) 


( 22 ) 


Now 


tan o * c/d 


(23) 


and 


c = (a 2 + b 2 ) 1 / 2 


(24) 


where 


a = b tan £ 


(25) 
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and 


d = b tan K 


(26) 


where a, b, c, and d are distances defined In figure 6, so 


tan a ■ b(tan 2 S + l)^ 2 /(b tan k) 


(27) 


or 


tan a « (tan 2 *- + l) 1/2 /tan k (28) 


Once the value of o can be determined, v and C can be evaluated In a man- 
ner similar to the general axl symmetric case where: 


cos v » 1/(1 + tan 2 cr + tan 2 v)^ 2 


(29) 


and 


COS C * (1 - COS 2 a - COS 2 v) 1/2 


(30) 


It Is possible to describe the Interception point and the surface normals 
at that point. Now the changes In probe volume geometry can be calculated. 


TOLERANCES 

The computer program requests user Input of four tolerances. There are 
two types of tolerances. The first type specifies how near parallel to a coor- 
dinate axis a direction vector need be to be considered parallel to that axis. 
This simplifies some calculations and can Improve accuracy by eliminating divi- 
sions by numbers near zero. The other type of tolerance Is used to specify 
the solution accuracy of the window surfaces during the calculations of the 
general smooth windows. 

All requested tolerances have default values. 

The first tolerance requested, 'TL1 1 , Is used during the calculation of 
the direction cosines of the nonrefracted laser beams. In later calculations, 
divisions by the direction cosines of the laser beams are required. The toler- 
ance ' TL1 ' specifies at what point the beam is considered to be parallel to a 
coordinate axis. This insures that division by a very small number during the 
calculation of the beam surface intercept point does not occur. 
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Tolerance ' TL2 ’ restricts the potential error generated in the solution 
of the general surface outward normal. Small errors in these values can make 
noticeable differences in the final solution. To limit the potential error, 
the slopes of the four points bracketting the proposed interception point are 
calculated. The changes in the slopes in the two coordinate directions are 
calculated and these changes are required to be less than ' TL2 . ' This effec- 
tively limits the maximum possible error in the surface normal of the inter- 
cept. If this convergence criteria is not met, a finer grid describing the 
window surface is generated and the calculation is repeated. 

The third tolerance, ' TL3 ' , restricts the possible error of the intercept 
point of the incoming light beam and the window surface. The Intercept point 
calculated using a linear interpolation between the points bracketting the 
intercept region is compared with the intercept point calculated using a cubic 
spline interpolation. If the difference is within the tolerance ' TL3 ' the pro- 
gram continues. If this criteria is not met, a finer grid describing the sur- 
face is generated and the intercept calculation is repeated. 

The last input tolerance, 1 TL4 ’ , is similar to the first since it speci- 
fies whether or not a direction cosine is near parallel to the X-axis or lies 
in the X to Y or X to Z planes. Those special cases allow simplified 
calculations and use of the simplified calculations insures that division by a 
very small number does not occur. 

All default tolerances result in sufficiently accurate solutions for the 
types of window in use at the Lewis Research Center. For comparison, note 
that the typical window at Lewis Research Center might have a change in slope 
of 2.3° per centimeter in a direction perpendicular to the incoming Laser Beam. 


EXAMPLES 

A number of test cases for the flat plate and cylindrical window were run 
and the results are available in reference 8. A comparison of output with an 
unpublished algorithm and also with hand calculations show an exact agreement 
with both the flat plate and simple cylinder test cases. 

Since the window shapes for the general axi symmetric and smooth general 
windows can only be approximated, ideal comparisons are difficult. In results 
not given here, the use of the general axisymmetric technique to calculate win- 
dow effects of a cylindrical window results in an identical solution with the 
cylindrical window solution and using the smooth general technique to calcu- 
late the effects of a flat plate window also results in an exact comparison 
with the flat plate solution. 

The general axisymmetric case was also run to calculate the window effects 
of an actual compressor window. Figure 7 shows the input geometry. The 
bisecting line is moved out of the X to Z and X to Y planes. This gen- 
eral axisymmetric case was also compared with two outputs run using the general 
smooth window algorithm. These two general smooth cases varied only with 
respect to the orientation of the window and lens geometry to the coordinate 
system. The first orientation run was similar to the general axisymmetric 
orientation shown in figure 7. The second orientation is shown in figure 8 
and simply results from rotating the problem geometry about the Y-axis. The 
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window geometry is given in Table 1(a) and the other relevant parameters are 
given in Table 1(b) . 

A comparison of outputs, provided in figures 9 to 11, showed very good 
agreement between the three cases. 

Figure 9 shows the uncrossing results, in fringes, versus S for the 
three cases. Maximum difference between cases is slightly over three fringes 
at S equal 90°. This is a relatively large number of fringes; however, it 
occurs in a geometry where the total uncrossing is a minimum. Over most of 
the range of comparison, the difference is less than two fringes. From a prac- 
tical standpoint, the agreement for all three cases is excellent. 

Figure 10 shows the change in beam crossing angled versus $ for the 
three test cases. Once again, the agreement is very good with the maximum dif- 
ference in the three configurations less than 0.017®. This maximum difference 
results in velocity errors of less than 0.3 meters per second. Note, however, 
that the maximum calculated velocity error due to the change in crossing angle 
from refraction effects is calculated to be about 1 meter per second for these 
test cases. 

Figure 11 shows the change in probe volume position for each of the three 
cases. The figure shows the relative movement in each coordinate direction 
and the total relative movement. Movement is defined as the change in probe 
volume position from the position it would occupy if there were no intervening 
window. Relative movement is that distance divided by the window thickness. 
Agreement is excellent for total movement and X direction movement with a 
maximum difference being less than 10 percent or, in absolute differences, 
0.00635 cm (window thickness of 0.3175 cm). The Y and Z direction move- 
ments show excellent agreement between the two test cases with the same win- 
dow/coordinate system orientation. However, when the window surface is 
rotated to an orientation that is not approximately in the Y to Z plane, 
the agreement is not as good. Nevertheless, for both directions, the total 
difference in movement between test cases is only about 0.00225x window thick- 
ness or about 7.15x1 0 - ^ cm for a 0.3175 cm thick window. This is equivalent 
to less than one third of the probe volume diameter. 

This set of test cases demonstrates two Important points. The first is 
to use the surface description that most accurately reflects the actual window 
geometry. If the window can be described using a general axi symmetric surface, 
this is the most desirable method of input. If the window can be specified as 
a flat plate or simple cylinder, these are the most desirable datasets. The 
second point is that, if the surface to be considered is best described as a 
smooth general surface, the orientation of the window surface normal should be 
kept as nearly parallel to the X-axis as possible. 


CONCLUSIONS 

A FORTRAN computer program has been developed to solve for the movement, 
beam uncrossing, and bisecting angle change of a probe volume of a laser fringe 
anemometer when it must traverse an Intervening region of different refractive 
index. 
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This program analytically solves for a flat plate window or a cylindrical 
window. It uses a cubic spline curve fit and an iterative procedure to solve 
for a general axi symmetric and a general surface window. 

The comparison of results for all window surface routines is good provided 
the surfaces are well behaved and well described in the region where the beams 
Intersect the window. 

The assumption that window refraction effects on the probe volume are neg- 
ligible cannot be supported in cases where complex window geometries are being 
used. 
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APPENDIX A 


INPUTS AND OPERATION OF THE PROGRAMS 

As currently configured, two separate computer codes are necessary for a 
complete solution. Both programs are Interactive and require user Inputs. 

Both programs use FORTRAN coding for the terminal Interactive operations and 
for manipulating the necessary data files. 

The first program sets up an Input data set to be run. This allows a per 
manent file of different input data sets to be maintained. The Input data 
file structure is shown in figure A1 . A glossary of terms is provided in 
Appendix B. 

Lines 4 through 6 are optional. Should nop2 be set to zero, the program 
will generate output for beam plane orientations every 10° from 0 to 170. In 
this Instance, no entries need to be made for lines 4 through 6 and the user 
will be Immediately prompted for line 7. 

Starting with line 9, the input varies depending on the nature of the win 
dow to be described. For a general window, the "upper" or outer surface varia 
bles are read in one at a time increasing in the ' Y ' or leading index first, 
as shown In figure 4. When finished with the window upper surface, continue 
with the lower or Inner surface. 

Inputs for the flat plate are done in the same manner except that only 
four points per surface (total of eight for the window) should be entered. 

For the cylindrical window, the only input past card 8 is the outer sur- 
face radius. 

For the general axi symmetric window, type in the number of Z to R 
points to describe the upper surface and the Z to R values. 

Figure A2 shows the coordinate system and variables used to specify the 
problem geometry. Angles should be input as degrees and distances can be spec 
ified In any consistent set of units. 

The second, or beam uncrossing program requires the file name containing 
the input data. 

The beam uncrossing program then queries the operator to Insure that the 
input parameters are correct. This also allows modification of all inputs 
other than the basic window geometry. 
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APPENDIX B 
INPUT DICTIONARY 


All Integer fields are 15 format and all real fields are F10 format. 

delta (l-nop2) an array of angles describing the angle between the 

plane the two incoming beams lie in and a plane contain- 
ing the lens center-z axis (axis of rotation), see fig- 
ure A2b, in degrees. 


D2XDY2 

eps 

iyi , iy2 

1 zl , iz2 

nl , n2 
nr 
nt 
nopl 

nop2 

rcl (1-nr) 


specifies the magnitude of the second derivative at the 
end points of the cubic spline fit when the window to 
be modeled is a general axi symmetric or general smooth 
window, (default*. 5) 

the bisecting angle of the beams in degrees 

number of points in the "y" direction to describe sur- 
faces 1, 2; 01 for cylindrical, 02 for flat plate, and 
10(max) for general surface 

number of points in the "z" direction to describe sur- 
faces 1, 2; 02 for flat plate, 10(max) for general 
surface 

index of refraction of the medium before and after trav- 
ersing the window; (default ■ 1.00029, air) 

number of Z to R coordinate pairs used to describe 
a general axi symmetric window 

index of refraction of the window; (default = 1.6, nomi- 
nal glass) 

flag to define the type of window to be Input; 0 - gen- 
eral surface, 1 - flat plate, 2 - cylindrical, and 3 - 
general axi symmetric 

defines whether the series of angles (delta) is 
defaulted or specified; 0 - defaulted, gives delta 
increments of 10° from 0 to 170, 1 to 18 indicates nop2 
number of delta angles to be input 

outer surface radius, general axisymmetric or cylindri- 
cal surface. 


rl 


distance from the center of the focusing lens to the 
axis of rotation (z axis) 


rpv 


distance of the probe volume (lens focal point assuming 
no intervening window) to the axis of rotation (z axis) 
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thk 


window normal thickness; not used for general surface 
input. 


thl angle between the x axis and a perpendicular line con- 

necting the z axis and the center of the focusing lens 

thpv angle between the x axis and a perpendicular line con- 

necting the z axis and the probe volume (focal point of 
the lens) 

til tolerance used in subroutine 'ppos.' Defines whether 

or not the plane in which the two Incoming laser beams 
lie is parallel to any coordinate plane. This speci- 
fies the nature of the calculation of the beam/focusi ng 
lens plane intercept calculation, (default is .0017rad) 

t 1 2 tolerance in subroutine 1 sf i n ' . Determines a success- 

ful solution to the general surface/laser beam inter- 
cept. This criteria restricts the maximum error 
between the actual surface normal direction and the cal- 
culated value. Default is set to cos (angle error)<.01. 

tl3 tolerance in subroutine ' sf 1 n ' . Determines a success- 

ful solution to the general surface/laser beam inter- 
cept. This is a limit between the difference in calcu- 
lated intercept points using the 'linint' (linear inter- 
polation) approximation of the window surface and the 
'splint 1 (spline interpolation) approximation of the 
window surface. It is assumed that when the difference 
is less than (tl3*thk), the surface description is accu- 
rate for an acceptable answer. Default is 0.0015. 

t 1 4 tolerance in subroutine ' sf i n ' . Determines if any sur- 

face direction vector is near parallel to the x-axis 
or liens in the X to Y or X to Z planes. Maxi- 
mum off parallel or out-of-plane error for the direc- 
tion cosine is tl4. Default is .0017. 

xl ,yl ,zl(l-1yl , 1 - 1 z 1 ) x, y, z coordinates describing the upper surface for 

the flat plate or general window case. 

x2,y2,z2(1-1y2,l-iz2) x, y, z coordinates describing the lower surface for 

the flat plate or general surface case 

Laser beam wavelength in Angstroms. Used for calculat- 
ing the uncrossing in fringes. (default=5145A) 

z coordinates (along the axis of symmetry) correspond- 
ing to rcl for the general axi symmetric case 

axial distance of focusing lens center 

axial distance of probe volume (focal point) 


wvl 

zcl (1-nr) 

zl 

zpv 
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APPENDIX C 


SYMBOLS 

f prove volume fringe crossing frequency 

n index of refraction 

N outward surface normal at point of laser beam/surface intersection 
V fluid (particle) velocity vector 
V, laser beam direction vector, 1 « 1 or 2 
c.fS.Y direction angle for incoming laser beam 

v,C,o direction angles for outward surface normal vector (fig. 5 and 6) 

X light wavelength 

< angle between a surface tangent line in the X to Z plane and the Z 

axis (fig. 5 and 6) 

5 angle between a surface tangent line in the X to Y plane and the 

Y axis (fig. 5 and 6) 

e beam crossing angle in probe volume 

$ angle between the outward surface normal and the Incoming laser beam 
Superscript 

' Indicates the region after the laser beam has crossed the surface 
boundary 

" indicates a coordinate system centered on the surface/ray Intercept 
point. Used in figures 4 and 6. 
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APPENDIX D 


SAMPLE INPUTS AND OUTPUTS 

Figures D1 through D4 show the sample Inputs for flat plate, simple cylin- 
der, general axisymmetric, and general smooth windows, respectively. In these 
figures, all user Inputs are underlined and an asterisk is found at the right 
hand margin of these lines. 

Notice that the data set up program allows the option to change any input 
data line. The line number of the line that is to be changed is typed in, 
then the entire line is retyped. 

Figure D5 shows a sample output from the main optics code. This output 
is self explanatory. 
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hello 


this is the optics data set-up program 

filename: ndat 1 

please input a dataset title. 

Flat Plate test case 1 

use a 4i5 field to type in nopl ,nop2,nun, iyl , izl , iy2, iz2. 
the surface options are specified by: 
nopl=0: general surface 
Isflat plate 
2: cylinder 

3: general axisymmetric 

the number of input "zt" angles is specified by nop2. default is 18. 
unit system is specified by nun (0=meters) (l=inches). 

nopl nop2 nun iyl izl iy2 iz2 

* * * * * * * 

1 0 1 2 2 2 2 

nopl= 1 nop2= 0 nun=*l iyl= 2 izl= 2 iy2= 2 iz2= 2 

ok, now type in rpv,zpv,thpv,rl,zl,thl 
use an 8fl0.5 field* 

* * * * * * * * 

8.8077 6.5334 3.6942 15. 3.4302 3.3399 

rpv= 8.80770 zpv= 6.53340 thetapv= 5.69420 

rl = 15.00000 zl = 3.43020 thetal = 3.33990 

ok, you have defaulted on the delta angles l 

type in the eps , thk,nl ,nt ,n2 ,d2xdy2,wvl in that order. 

default for d2xdy2 is .5 and for wvl it is 5145A. 

* * * * * * * 

-> 11.49 .125 1.00029 1.52 1.00029 .5 5145 . 

epsilon= 11.490 thickness^ 0.1250 ni= 1.00029 nt= 1.52000 n2= 1.00029 

splint d2xdy2= 0.50000 wvl = 05145. 

type in the tolerances til thru tl4,use a 4fl0 field. 

* * * * * * * 


for a flat plate, you must describe each surface 
with four x,y,z, points. 

now you must type in zl ,yl ,xl , then z2,y2,x2 in the 
order of increasing iy then iz; use an 8fl0 field. 

this general surface has 8 input points 

* . * * * * * *. 

CK 0. 2^ 

* — — * * * * *. 

4 K 2^ 

* * - * — * * *— — 

* 1 . 0 ^ h. 

* * * * * * *. 

-> 1. U 

* * * * - — *- * *. 

0. Ch 1.875 

* * * * * * *. 

* 0. 1_. 1.875 

* * * * . * * *• 

U 0^ 1.875 

* * * * - — * -* 

-» U 1. 1.875 


FIGURE Dl. - FLAT PLATE INPUT DATASET EXAMPLE CASE. ARROWS INDICATE 
USER INPUT. 


itA 



1 0 

12 2 2 

2 




file line 

for the next line is 

3. 




8.80770 

6.53540 5.69420 ' 

15. 00000 

3.43020 

3.33990 

0.00000 

file line 

for the next line is 

4. 




11.49000 

0.12500 1.00029 

1.52000 

1.00029 

0.50000 

0.05145 

file line 

for the next line is 

5. 




0.00000 

0.00000 0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

6. 




0.00000 

0.00000 2.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

7. 




0.00000 

1.00000 2.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

8. 




1.00000 

0.00000 2.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

9. 




1.00000 

1.00000 2.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

10. 




0.00000 

0.00000 1.87500 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

11. 




0.00000 

1.00000 1.87500 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

12. 




1.00000 

1.00000 1.87500 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

13. 




1.00000 

1.00000 1.87500 

0.00000 

0.00000 

0.00000 

0.00000 

if you wish to change an input 

line; type 

in the ! 

Line . 



use an i3 field, else default. 
* 


FIGURE Dl. - CONCLUDED. 


I SB 


0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 



hello 


this is the optics data set-up program 

filename: ndat2 

please input a dataset title. 

4 Simple Cylinder test case 1 

use a 4i5 field to type in nopl ,nop2,nun,iyl ,izl ,iy2,iz2. 
the surface options are specified by: 
nopl=0: general surface 
l:flat plate 
2: cylinder 

3: general axisymmetric 

the number of input "zt" angles is specified by nop2. default is 18. 
unit system is specified by nun (0=meters) (l*inches). 

nopl nop2 nun iyl izl iy2 iz2 

* * * * * * * 

2 0 1 1 . 

nopl* 2 nop2= 0 nun=l iyl= 1 izl* 0 iy2= 0 iz2= 0 

ok, now type in rpv,zpv,thpv,rl,zl,thl 
use an 8fl0.5 field. 

* * it * it * * it 

-> 8.8077 6.5354 5.6942 15. 3.4302 3.3399 

rpv* 8.80770 zpv* 6.53540 thetapv* 5.69420 

rl = 15.00000 zl * 3.43020 thetal * 3.33990 

ok, you have defaulted on the delta angles! 

type in the eps, thk,nl ,nt f n2,d2xdy2,wvl in that order. 

default for d2xdy2 is .5 and for wvl it is 5145A. 

it * it it * it * 

11.49 .125 1.00029 1.52 1.00029 

epsilon* 11.490 thickness* 0.1250 nl= 100029 nt* 1.52000 n2= 1.00029 

splint d2xdy2= 0.00000 wvl* 0.00000 

type in the tolerances til thru tl4,use a 4fl0 field. 

* it it * it it * 


you have selected a cylindrical window. 

you must type in the radius of the outer cylinder 

surface, use an flO field. 

★ * * it * * it. 

12 

outer cylinder radius is: 10.00000 


**** this is the read file data **** 
file line for the next line is 3. 


8.80770 

6.53540 

5.69420 

15.00000 

3.43020 

3.33990 

0.00000 

0.00000 

file line 

for the next line is 

4. 





11.49000 

0.12500 

1.00029 

1.52000 

1.00029 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

5. 





0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

6 . 





10.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 


if you wish to change an input line; type in the line, 
use an i3 field, else default. 

* 


FIGURE D2. - CYLINDRICAL INPUT DATASET EXAMPLE CASE. ARROWS 
INDICATE USER INPUT. 


JfrC. 



hello 


this is the optics data set-up program 
filename: ndat!3 

use a 4i5 field to type in nopl,nop2,nun,iyl,izl,iy2,iz2. 
the surface options are specified by: 
nopl=0: general surface 
l:flat plate 
2: cylinder 

3: general axisymmetric 

the number of input "at” angles is specified by nop2. default is 18. 
unit system is specified by nun (0«meters) (l*inches). 

nopl nop2 nun iyl izl iy2 iz2 

* * * * rt * * 

-> 3 0 1 * 10 10 10 10 

nopl= 3 nop2= 0 nun=l iyl= 10 izl= 10 iy2= 10 iz2= 10 

ok, now type in rpv,zpv, thpv,rl,zl,thl 
use an 8fl0.5 field. 

it * * * it it it * 

-> 8.8077 6.5354 5.6942 13. 3.4302 3.3399 

rpv= 8.80770 zpv* 6.53540 thetapv= 5.69420 

rl * 15.00000 zl = 3.43020 thetal = 3.33990 

ok, you have defaulted on the delta angles! 
type in the eps,thk,nl,nt,n2,d2xdy2,wvl. 

★ it it * it 

-> 11.49 .125 1.00029 1.52 1.00029 .5 

epsilon* 11.490 thickness* 0.1250 nl* 100029 nt= 1.52000 wvl* 0.00000 
type in the tolerances til thru tl4,use a 4fl0 field. 

* it * it * 


you have selected a general axisymmetric surface 

use an i3 field to specify the number of outer surface points (10 max) 

— it 

-» 10 

use a 2fl0 field to read in z then r 

* it * it * 



3.7 

10.23 


4.2 

10.225 


5.1 

10.217 


6. 

10.205 


6.727 

10.125 


7.3 

9.995 


7.7 

9.915 


8.2 

9.842 


8.7 

9.791 


9.2 

9.755 


FIGURE D3. - GENERAL AXISYMMETRIC INPUT DATA SET EXAMPLE CASE. 
ARROWS INDICATE USER INPUT. 
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**** this is the read file data **** 


3 0 

1 10 10 10 

10 





file line 

for the next line is 

2. 





8.80770 

6.53540 5.69420 

15.00000 

3.43020 

3.33990 

0.00000 

0.00000 

file line 

for the next line is 

3. 





11.49000 

0.12500 1.00029 

1.52000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

4. 





0.00000 

0.00000 0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

5. 





3.70000 

10.23000 0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

6. 





4.20000 

10.22500 0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

7. 





5.10000 

10.21700 0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

8. 





6.00000 

10.20500 0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

9. 





6.72700 

10.12500 0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

10. 





7.30000 

9.99500 0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

11. 





7.70000 

9.91500 0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

12. 





8.20000 

9.84200 0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

13. 





8.70000 

9.79100 0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

14. 





9.20000 

9.75500 0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 


if you wish to change an input line; type in the line, 
use an i3 field, else default, 

— * 

FORTRAN STOP 


FIGURE D3. - CONCLUDED. 




hello 


4 $ run opt in 


this is the optics data set-up program • 

4 filename: ndat6 

please input a dataset title. 

4 NDAT6 GENERAL SURFACE INPUT 

use a 4i5 field to type in nopl,nop2,nun,iyl,izl,iy2, iz2. 
the surface options are specified by: 
nopl=0: general surface 
l:flat plate 
2: cylinder 

3: general axisymmetric 

the number of input "zt" angles is specified by nop2. default is 18. 
unit system is specified by nun (0=meters) (l=inches). 

nopl nop2 nun iyl izl iy2 iz2 

* * * * * * * 

4 0 0 1 5 5 5 5 

nopl= 0 nop2= 0 nun=l iyl= 5 izl= 5 iy2= 5 iz2= 5 

ok, now type in rpv,zpv, thpv,rl ,zl , thl 
use an 8fl0.5 field. 

* * * * * * * *. 

4 8.8077 6.5354 5.6942 15. 3.4302 3.3399 

rpv* 8.80770 zpv= 6.53540 thetapv= 5.69420 

rl = 15.00000 zl = 3.43020 thetal = 3.33990 

ok, you have defaulted on the delta angles! 
type in the eps , thk,nl ,nt ,n2,d2xdy2,wvl in that order, 
default for d2xdy2 is .5 and for wvl it is 5145A. 


* - * * — * 

4 11.49 .125 1.00029 1.52 

* - 

1.00029 1. 



epsilon= 11.490 thickness= 0.1250 nl= 
splint d2xdy2= 1.00000 wvl* 0.00000 
type in the tolerances til thru tl4,use 
* * * * 

100029 nt= 

a 4fl0 field. 
* * 

1.52000 n2= 

1.00029 

4 





now you must type in zl,yl,xl,then z2,y2,x2 in the 
order of increasing iy then iz, use an 8fl0 field. 


this general surface has 50 input points 


4 

5.624959950.5346813810.2023287 

* 

* 

* 


4 

5.624959950.7126551310.1914434 

* * * *.._WT 

* 


* 

4 

5.624959950.8904118510.177454 





FIGURE D4. - GENERAL SMOOTH INPUT DATASET EXAMPLE CASE. ARROWS 
INDICATE USER INPUT. 


I TF 



6 . 144859790. 7025146510.0464277 



* 

* 

* _ 

6.144859790.8777419910.0326366 

* _ * * _ * 

•ft — 

* 

..... * .... 

6.144859791.0527019510.01579 
* * * * 

* 

* 

* 

6.144859791.227341299.99589252 


0 0 

15 5 5 

5 





file line 

for the next line is 

3. 





8.80770 

6.53540 5.69420 

L5. 00000 

3.43020 

3.33990 

0.00000 

0.00000 

file line 

for the next line is 

4. 





11.49000 

0.12500 1.00029 

1.52000 

1.00029 

1.00000 

0.00000 

0.00000 

file line 

for the next line is 

5. 





0.00000 

0.00000 0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

6. 





5.62496 

0.53468 10.20233 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

7. 





5.62496 

0.71266 10.19144 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

8. 





5.62496 

0.89041 10.17745 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

9. 





5.62496 

1.06790 10.16036 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

10. 





5.62496 

1.24506 10.14018 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

11. 





5.73771 

0.53458 10.20044 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

12. 





5.73771 

0.71252 10.18956 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

39. 





5.80099 

1.05447 10.03261 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

40. 





5.80099 

1.22940 10.01268 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

41. 





5.91475 

0.52775 10.07014 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

42. 





5.91475 

0.70342 10.05940 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

43. 





5.91475 

0.87888 10.04559 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

44. 





5.91475 

1.05406 10.02872 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

45. 





5.91475 

1.22893 10.00880 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

46. 





6.02933 

0.52746 10.06460 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

47. 





6.02933 

0.70303 10.05386 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

48. 





6.02933 

0.87839 10.04006 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

49. 





6.02933 

1.05348 10.02320 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

50. 





6.02933 

1.22825 10.00329 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

51. 





6.14486 

0.52707 10.05716 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

52. 





6.14486 

0.70251 10.04643 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

53. 





6.14486 

0.87774 10.03264 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

54. 





6.14486 

1.05270 10.01579 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

file line 

for the next line is 

55. 





6.14486 

1.22734 9.99589 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 


if you wish to change an input line; type in the line, 
use an i3 field, else default. 

— * 

FORTRAN STOP 


FIGURE D4. - CONCLUDED. 
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******************************************************************** 

THIS IS THE OUTPUT FOR THE OPTICS UNCROSSING PROGRAM 

******************************************************************** 

********** INPUT data ********** 


Original file name was NDAT13 

Window surfaces are "general axisymraetric surf aces", (nol=3) 

The input data is provided in English units (inches ,nun=l ) 

The desired 18 input beam plane orientation angles, in degrees, are: 
0.000 10.000 20.000 30.000 40.000 50.000 60.000 70.000 80.000 
90.000 100.000 110.000 120.000 120.000 130.000 140.000 150.000 170.000 


The probe volume and lens geometry are: 

z 

axis 
6.5354 
3.4302 

The Xi (out-of -plane ) angle is -26.56548 


pv: 

lens: 


radial axial angular x y 

position position position axis axis 

8.8077 6.5354 5.6942 8.7642 0.8739 

15.0000 3.4302 3.3399 14.9745 0.8739 


Indices of refraction are: lens-side medium nl= 

window medium nt= 
pv-side medium n2= 


1.00029 

1.52000 

1.00029 


Original beam crossing angle is 11.49 deg 


Window thickness is 0.1250 

Laser wavelength is 5145A giving a fringe spacing of 0.101181E-03 length/fringe 
Input tolerances have been defaulted. 

The General Axisymmetric upper and lower surface points are: 


axial 

upper 

lower 

losition 

radius 

radius 

3.700 

10.230 

10.105 

4.200 

10.225 

10.100 

5.100 

10.217 

10.092 

6.000 

10.205 

10.080 

6.727 

10.125 

9.998 

7.300 

9.995 

9.867 

7.700 

9.915 

9.788 

8.200 

9.842 

9.716 

8.700 

9.791 

9.666 

9.220 

9.755 

9.630 


FIGURE D5. - SAMPLE OUTPUT 
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******************************************************************* 

********** OUTPUT DATA ********** 

probe volume movement 

definitions .’delta .’beam plane orientation angle 

xnew,ynew,znew: the new probe volume position 

xrel »yrel ,zrel : the change in probe volume position /window thickness 
i.e. (Xpv-Xorig)/thk 

reltot :the vector sum of the change in probe volume 

position/window thickness 


delta 

xnew 

ynew 

znew 

xrel 

yrel 

zrel 

reltotal 

(deg) 

(inches) 

(inches) 

(inches) 





0.00 

8.71146 

0.87040 

6.54039 

0.42226 

0.02786 

-0.03989 

0.42505 

10.00 

8.71176 

0.87039 

6.54025 

0.41982 

0.02796 

-0.03879 

0.42253 

20.00 

8.71316 

0.87039 

6.53960 

0.40860 

0.02801 

-0.03358 

0.41093 

30.00 

8.71548 

0.87038 

6.53852 

0.39004 

0.02805 

-0.02497 

0.39184 

40.00 

8.71838 

0.87038 

6.53717 

0.36684 

0.02806 

-0.01420 

0.36819 

50.00 

8.72151 

0.87038 

6.53572 

0.34181 

0.02804 

-0.00258 

0.34297 

60.00 

8.72448 

0.87039 

6.53435 

0.31811 

0.02800 

0.00842 

0.31945 

70.00 

8.72694 

0.87040 

6.53320 

0.29838 

0.02794 

0.01758 

0.30020 

80.00 

8.72864 

0.87040 

6.53241 

0.28480 

0.02788 

0.02389 

0.28715 

90.00 

8.72939 

0.87041 

6.53206 

0.27879 

0.02780 

0.02670 

0.28144 

100.00 

8.72912 

0.87042 

6.53218 

0.28093 

0.02774 

0.02573 

0.28347 

110.00 

8.72786 

0.87043 

6.53276 

0.29102 

0.02769 

0.02108 

0.29309 

120.00 

8.72574 

0.87043 

6.53375 

0.30801 

0.02765 

0.01323 

0.30954 

130.00 

8.72297 

0.87043 

6.53503 

0.33014 

0.02764 

0.00298 

0.33131 

140.00 

8.71987 

0.87043 

6.53647 

0.35494 

0.02765 

-0.00853 

0.35612 

150.00 

8.71681 

0.87043 

6.53789 

0.37944 

0.02769 

-0.01991 

0.38097 

160.00 

8.71417 

0.87042 

6.53912 

0.40053 

0.02774 

-0.02973 

0.40259 

170.00 

8.71230 

0.87041 

6.53999 

0.41550 

0.02781 

-0.03672 

0.41804 


FIGURE D5. - CONTINUED. 
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Page 19 

******************************************************************* 


probe volume uncrossing and crossing angle changes 


definitions sancross: crossing angle 

delang ; change in crossing angle/original angle in percent 
uncross: total uncrossing in inches or meters 
relun : uncrossing/window thickness 
uncrosf : uncrossing in fringes 


delta 

ancross 

delang 

uncross 

relun 

uncrosf 

(deg) 

(deg) 

(percent) 

(inches) 


(fringes) 

0.00 

11.44312 

0.404E+00 

0. 168930E-03 

0. 135144E-02 

1.670 

10.00 

11.44699 

0. 370E+00 

-.522350E-03 

-.417880E-02 

-5.163 

20.00 

11.45611 

0. 291E+00 

-.114614E-02 

-. 916913E-02 

-11.328 

30.00 

11.46944 

0. 175E+00 

-. 162182E-02 

-. 129746E-01 

-16.029 

40.00 

11.48513 

0.384E-01 

-. 189497E-02 

151598E-01 

-18.729 

50.00 

11.50134 

-0, 103E+00 

-.193423E-02 

154739E-01 

-19.117 

60.00 

11.51596 

-0. 230E+00 

-.174236E-02 

-. 139389E-01 

-17.220 

70.00 

11.52743 

-0.330E+00 

-. 134613E-02 

-. 107691E-01 

-13.304 

80.00 

11.53436 

-0 . 390E+00 

-.796640E-03 

-.637312E-02 

-7.873 

90.00 

11.53601 

-0.404E+00 

-. 158393E-03 

-. 126714E-02 

-1.565 

100.00 

11,53226 

-0.372E+00 

0.498078E-03 

0. 398462E-02 

4.923 

110.00 

11.52354 

-0 . 296E+00 

0 , 109834E-02 

0.878668E-02 

10.855 

120.00 

11.51076 

-0.185E+00 

0. 157398E-02 

0. 125919E-01 

15.556 

130.00 

11.46554 

-0.522E-01 

0. 186899E-02 

0 . 149519E-01 

18.472 

140.00 

11,47945 

0.878E-01 

0.194407E-02 

0. 155525E-01 

19.214 

150.00 

11,46446 

0. 218E+00 

0. 178272E-02 

0. 142617E-01 

17.619 

160.00 

11.45249 

0. 322E+00 

0. 139910E-02 

0. 111928E-01 

13.828 

170.00 

11,44504 

0. 387E+00 

0.838048E-03 

0.670438E-02 

8.283 


FIGURE D5. - CONCLUDED. 
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APPENDIX E 


WARNING MESSAGES 

A number of warning messages may be printed during the course of a calcu- 
lation. They are, for the most part, only advisory In nature with the calcula- 
tion continuing. However, they may Indicate an error In the input data. The 
possible warning messages are: 

1. "WARNING: FROM RTRACE , INTERNAL REFLECTION OCCURS IN THIS CONFIGURATION." 

It Is possible for an optics window configuration to occur where the beam leav- 
ing the window intercepts the window surface at an angle greater than Bragg's 
angle. Should this occur, an Internal reflection of the beam will occur. The 
calculation stops and output from this configuration should be considered 
Incorrect. This situation should only occur with general smooth and general 
axi symmetric windows. 

2. "WARNING: FROM LINT, LINES DO NOT CROSS IN SPACE." 

This warning is printed out when calculations indicate that the two incoming 
beams lie in the same plane but intercept point calculations do not indicate 
the beams cross. This message indicates an input data error. 

3. "WARNING: FROM CYLINT, LINE DOES NOT INTERCEPT CYLINDER WALL." 

This message indicates an input geometry where the beams do not intercept the 
cylinder window outer surface. This message indicates an input data error. 
Calculation terminates. 

4. "WARNING: FROM CYLINT, THE CALCULATIONS ON Z POSITION DO NOT AGREE." 

This message indicates an input data error or severe geometry constraints such 
as an intercept near the surface tangency point or an excessive 4>. Calcula- 
tion continues. 

5. "WARNING: FROM CYLINT, CALCULATED INTERCEPT X AND Y DOES NOT LIE ON R. 

ALTERNATE CALCULATION STARTED." 

This message indicates severe geometry constraints where the surface intercept 
point lies near the surface tangency point. The alternate calculation should 
provide greater accuracy for this geometry. Calculation continues. 

6. "WARNING: FROM AXI , NO CONVERGENCE IN 'AXI'." 

The calculation of the intercept point of the general axisymmetric surface and 
an incoming ray is an iterative procedure. Should the difference between the 
calculated line intercept and the surface intercept point be greater than TL1 
after 50 iterations, this message is printed. Calculation continues; however, 
the output data may be in error. This message indicates an input data error 
or severe surface curvature problems. 
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7. "WARNING: FROM SURLNT, SURLNT FAILS." 

This message indicates that the intercept of the incoming light ray with the 
linear approximation to the general surface fails to converge within 76 itera- 
tions. Calculation continues; however, a possible problem may exist with the 
input surface data. 

8. "WARNING: FROM SFIN, CHANGES IN THE SURFACE OR THE SURFACE NORMAL ARE 

RAPID, SFIN DOES NOT CONVERGE. SOLUTION ACCURACY IS NOT GUARANTEED." 

This message indicates a severe window geometry resulting in a solution which 
may be in error. Calculation continues. 

9. "WARNING: FROM FINEE, MORE THAN ONE DIRECTION COSINE IS EQUAL TO ONE." 

This message should not occur and indicates an input data error that is proba- 
bly in the lens-PV geometry. Calculation terminates. 

10. "WARNING: FROM FINEE, THE BEAMS ARE PARALLEL AND DO NOT CROSS." 

This message indicates that the window geometry has resulted in two parallel 
beams. There will be no probe volume and no minimum crossing distance. Calcu- 
lation terminates. 

11. "WARNING: FROM FINEE, DMIN-' DMN* 1 

The minimum crossing distance is calculated in two ways. Should the results 
be different, this advisory only message is printed. A large difference in 
values would indicate an input data error. 
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TABLE I 


(a) General axi symmetri c window 
geometry 


Axial Position, 
cm 


0 . 

9.398 

10.668 

12.954 

15.24 

16.386 

17.0866 

18.288 

18.542 

19.558 

19.947 

20.828 

22.098 

23.368 

24.638 

25.908 

28.245 


Outer wi ndow 
radial position, 
cm 


25.9842 

25.9842 

25.9779 

25.9715 

25.921 

25.8159 

25.7228 

25.4716 

25.3949 

25.1892 

25.1367 

25.001 

24.8721 

24.7783 

24.7225 

24.7015 

24.7015 


(b) Lens probe volume configuration, general 
axi symmetric case 


Out of plane angle (Xi), deg 26.56548 

Probe volume position: 

radius, cm 22.3726 

axial position, cm 16.5999 

angular position, deg 5.6942 

Lens position: 

radial position, cm 38.1 

axial position, cm 87.12708 

angular position, deg 3.3399 

Original (measured) beam crossing angle, cm . . . 11.49 

Window normal thickness, cm 0.3175 

Window index of refraction 1.52 

Laser wavelength 5145A 


TABLE II. - GENERAL SURFACE TEST CASES 

[Comparison between calculated and measured 
test case probe volume movement and 
crossing angle change.] 


General axi symmetric test case 

Position change 

Cal culated 

Measured 

X 

0.032255 

0.0336 

Y 

.002665 

.0022 

2 

.00907 

.00685 

Crossing angle 
change 

.0068 

.001276 


General smooth test case 

Position change 

Cal cul ated 

Measured 

X( cm) 

0.04464 

0.0336 

Y ( cm) 

.00354 

.0022 

Z{ cm) 

.00011 

.000685 

Crossing angle 
change 

.02054 

.001276 








FIGURE 1. - FLAT PLATE GEOMETRY AND NOMENCLATURE. 
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FIGURE 3. - GENERAL AX I SYMMETRIC SURFACE AND IT'S COORDINATE SYSTEMS. 



FIGURE 2. - SIMPLE CYLINDER GEOMETRY AND NOMENCLATURE. 


r OUTER OR 



FIGURE 4. - GENERAL WINDOW INPUTS. 




(B) SURFACE NORMAL GEOMETRY AT POINT "P". 

FIGURE 5. - GENERAL AX I SYMMETRIC SURFACE GEOMETRY. 
NOTE: SUPERSCRIPTS INDICATE TRANSLATED AXES. 




(A) SURFACE NORMALS AND SLOPES AT SURFACE (POINT "P w ). 



(B) SURFACE NORMAL GEOMETRY AT POINT "P". 

FIGURE 6. - GENERAL SURFACE GEOMETRY. NOTE: SUPERSCRIPT IN- 
DICATES TRANSLATED AXES. 


LENS 

CENTER 

POSITION 



AXIAL 

POSITION 


(A) TEST CASE GEOMETRY, R-Z (X-Z) PLANE. 



(B) TEST CASE GEOMETRY, R-0 (X,Y) PLANE. 

FIGURE 7. - GENERAL AX I SYMMETRIC TEST CASE GEOMETRIES. 



FIGURE 8. - GENERAL SMOOTH TEST CASE GEOMETRY NUMBER 2 


BEAM UNCROSSING (FRINGES) 



FI6URE 9. - BEAM UNCROSSING IN FRINGES VERSUS 6, GENERAL AXISYMMETRIC AND GENERAL SMOOTH 
TEST CASES. 




FIGURE 10. - CHANGE IN BEAM CROSSING ANGLE VERSUS 6, GENERAL AXISYMMETRIC AND GENERAL SMOOTH 
TEST CASES. 
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O NDAT 6: GENERAL SURFACE, TILTED LOOK ANGLE 
□ NDAT 7: GENERAL SURFACE, TILTED SURFACE 
O NDAT 8: GENERAL AXISYMMETRIC SURFACE 
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(A) TOTAL RELATIVE PROBE VOLUME MOVEMENT VERSUS 
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(A) WINDOW GEOMETRY AND VARIABLES. 


(C) GENERAL AXISYMMETRIC. 
FIGURE A1 . - INPUT DATA FORMAT. 



(B) LENS/PROBE VOLUME GEOMETRY. 
FIGURE A2. - SYSTEMS GEOMETRY. 
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